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Abstract 



Using the D0 detector, we have observed events produced in pp coUisions that 
contain W or Z bosons in conjunction with very Uttle energy deposition ("rapidity 
gaps") in large forward regions of the detector. The fraction of W boson events 
with a rapidity gap (a signature for diffraction) is 0.89 iLq;];! %, and the probabihty 
that the non-diffractive background fluctuated to yield the observed diffractive 
signal is 3 X 10~^^, corresponding to a significance of 7.5a. The Z boson sample 
has a gap fraction of 1.44±o;52 with a significance of 4.4cr. The diffractive events 
have very similar properties to the more common non-diffractive component. 

Inelastic diffractive collisions are responsible for about 15% of the pp total cross section, and 
have been described by Regge theory through the exchange of a pomeron [1]. Such events are 
characterized by a proton (or antiproton) carrying away most of the beam momentum, and by 
the absence of significant hadronic particle activity over a large region of pseudorapidity (r] = 
— ln[tan(|)], where 9 is the polar angle relative to the beam). This empty region is called a rapidity 
gap and can be used as an experimental signature for diffraction. Ingelman and Schlein proposed 
the possibility of a partonic structure for the pomeron, which would lead to hard scattering in 
diffractive events [2]. This so-called "hard diffraction" was first observed by the UA8 experiment [3] 
at the CERN SppS collider in the form of jet events with an associated tagged proton. 

Initial rapidity-gap-based analyses of diffractive jet [4-6], 6-quark [7], and J /"if [8] production 
are qualitatively consistent with a predominantly hard gluonic pomeron, but the production cross 
sections observed at the Fermilab Tevatron are far lower than predictions based on data from 
the DESY ep collider HERA [4,9]. Diffractive jet results from the CDF collaboration using an 
antiproton tag [10] confirm the normalization discrepancy between Tevatron (^/s = 1800 GeV) 
and HERA data, while recent D0 rapidity-gap-based diffractive jet results at = 1800 GeV 
and y/s — 630 GeV [1 1] show that a simple normalization difference cannot accommodate the 
Tevatron data (and imply that a significant soft gluon component is needed to "save" the Ingelman- 
Schlein model). A unified picture within this framework requires a detailed understanding of 
gap survival probability, which includes effects from multiple parton scattering and extra gluon 
emission associated with the hard sub-process [12]. The soft color interaction (SCI) model [13], 
which hypothesizes that non-perturbative gluon emissions can create rapidity gaps, provides an 
alternative description of diffraction without invoking pomeron dynamics, and predicts diffractive 
rates similar to those observed. 

Bruni and Ingelman [14] proposed that a search for diffractive production of W and Z bosons 
would provide important information on diffractive structure, due to their sensitivity to quark 
sub-structure. They predicted that a pomeron composed primarily of quarks would lead to more 
than 15% of W and Z bosons being diffractively produced. The SCI model, on the other hand, 
predicts a diffractive fraction of about 1% [15]. 

The CDF collaboration observed a 3.8 standard deviation {a) signal for diffractive W boson 
production, extracting the signal using the asymmetry of both lepton charge and position relative 
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to the region of the rapidity gap, and obtained a diffractive to non-diffractive production ratio of 
(1.15 ± 0.55)% [16]. In this Letter, we present a definitive observation of diffractively produced 
W and Z bosons. We present characteristics of diffractive W bosons, and measurements of the 
fraction of W and Z boson events that contain forward rapidity gaps. In addition, we provide the 
ratio of diffractive W and Z cross sections, and the fraction of the initial momentum carried away 
by the scattered proton in the collision. 

In the D0 detector [17], electrons are measured and missing transverse energy {^x) deter- 
mined using the uranium/hquid- argon calorimeters, with electromagnetic coverage to |?7| = 4.1 
and coverage for hadrons to |?7| = 5.2. Electron identification, described in more detail below, 
requires a central or forward drift chamber track to match the location of the associated electro- 
magnetic cluster. For the period during which the data were collected, the D0 detector had no 
magnetic field within the central tracking volume, consequently electrons and positrons could not 
be differentiated and are both referred to as electrons. 

To identify rapidity gaps, we count the number of tiles containing a signal in the L0 forward 
scintillator arrays (tzlo) and the number of calorimeter towers {At] x = 0.1 x 0.1) with signals 
above threshold (tical)- The L0 arrays provide partial coverage in the region 2.3 < \r]\ < 4.3. 
A portion of the two forward calorimeters (3.0 < {rjl < 5.2) is used to measure the calorimeter 
multiplicity, with a particle tagged by the deposition of more than 150 (500) MeV of energy in 
an electromagnetic (hadronic) tower. These thresholds are set to minimize noise from radioactive 
decays in the uranium, while maximizing sensitivity to energetic particles [18]. 

For this analysis, we search for the presence of rapidity gaps in inclusive samples of — > ei/ 
events and Z —>■ e^e~ events, based on data with an integrated luminosity of approximately 
85 pb~^ accumulated during the 1994-1995 collider run (Run lb). The D0 collaboration has 
extensively studied W and Z boson production in the electron channel [19,20]. The requirements 
for the event selection in this analysis are nearly identical to those of Ref. [20], with two notable 
exceptions detailed below. The data were obtained using a single hardware trigger that required 
at least one electromagnetic (EM) object with transverse energy {Et) greater than 15 GeV, with 
more than 85% of its energy deposited in the EM section of the calorimeter (EM fraction). At 
the software trigger level, the EM cluster is required to satisfy isolation, shower-shape, and EM 
fraction criteria consistent with the presence of an electron. For the W boson sample, we require 
this candidate electron to have an Et > 20 GeV, and additionally require ^j, > 15 GeV for the 
neutrino, while for the Z boson sample, we require two electron candidates with Et > 16 GeV. 

The first significant difference between the data samples in this analysis and those of Ref. [20] 
is that we are unable to include events from the first portion of Run lb, during which a coincidence 
(in the L0 detector) between the remnants of the proton and antiproton was required, effectively 
vetoing single-diffractive production. Restricting this analysis to the part of the data collected 
without this condition reduces the sample by 30%. The only other major difference is that this 
analysis requires the removal of events with more than one proton-antiproton interaction in the 
same bunch crossing. This "single interaction" requirement is necessary for rapidity-gap-based 
diffractive studies, because the presence of additional events obscures the rapidity-gap signature. 
About 70% of the remaining data sample is discarded as a result of this requirement, which makes 
use primarily of timing information in the luminosity counters and the number of vertices found 
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in the central tracker to reject multiple interaction events. 



Variable 


Comment 


Events 


Trigger 


electron + 


119,890 


No L0 requirement in trigger 




84,310 


Main ring cuts 




63,978 


Single interaction 




17,870 


One electron in fiducial region 


\r)\ < 1.1 or 1.5 < \r]\ < 2.5 


17,626 


Et of electron 


> 25 GeV 


15,203 


Electron quality 


isolation, shape, EM fraction 


13,770 




> 25 GeV 


12,622 


Total W ^ eu sample 




12,622 


Central e sample 


\ri\ < 1.1 


8,724 


Forward e sample 


1.5 < \r]\ < 2.5 


3,898 



TABLE I. Central and forward W boson event selection criteria. 



The other analysis cuts are all standard criteria employed in D0 electron analyses. In addition 
to a 25 GeV threshold on the event and the electron Et, and further selection based on the 
electron quality, events that occurred during the injection of proton bunches in the Main Ring 
accelerator are rejected (these often produced significant energy deposition in the D0 calorime- 
ter) [19]. The final data samples consist of 811 Z boson candidate events, and 12,622 W boson 
candidate events, of which 8,724 have a central electron < 1.1) and 3,898 have a forward 
electron (1.5 < |?7| < 2.5). A summary of the event selections is given in Tables I and II. 

Figure 1 shows two views of nL0 versus ncAL for the combined central and forward W boson 
sample. The multiplicity in the forward rj interval with the lower ^cal multiplicity (for some 
events this interval is at +77 and for others —77) is plotted for Fig. 1(a) and (b) for the full 
range of multiplicity and for the region of low multiplicity (ncAL < 20, nL0 < 10), respectively. 
The distributions peak at zero multiplicity (ncAL = ^l0 = 0), in qualitative agreement with 
expectations for a diffractive component in the data. Figure 2 shows this scatter plot separately 
for the (a) central and (b) forward W boson samples, and for the (c) Z boson sample. All samples 
show clear evidence for a diffractive component at low multiplicity. 

We now compare characteristics of the diffractive W boson candidates to the non-diffractive 
events to verify that these are typical W bosons, except for the presence of a rapidity gap. Fig- 
ures 3(a), (c), and (e) show the electron Et, $t-i transverse mass (Mr), respectively, for 
standard W boson events (ncAL > 1), while Figs. 3(b), (d), and (f) show the corresponding quan- 
tities for diffractive candidate events (nL0 = ^cal = 0). Although the statistics in the diffractive 
sample are limited, the distributions in all three variables are very similar. The mean values for 
these quantities for the non-diffractive and diffractive samples, respectively, are in excellent agree- 
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Variable 


Comment 


Events 


Trigger 


two electrons 


13,912 


No L0 requirement in trigger 




10,023 


Main ring cuts 




8,751 


Single interaction 




2,381 


Two electrons in fiducial region 


\ri\ < 1.1 or 1.5 < |r/| < 2.5 


1,617 


Et of electrons 


> 25 GeV 


1,046 


Electron quality 


isolation, shape, EM fraction 


893 


Invariant mass window 


76 < Mee < 106 GeV/c^ 


811 


Total Z ^ e+e~ sample 




811 



TABLE II. Z boson event selection criteria. 




(a) <^cfrJ- (b) nc^^- 



FIG. 1. The multiplicity in forward tiles (nL0) and in calorimeter towers (ncAL) for the forward 
region with the lower multiplicity, for the combined central and forward electron W boson samples: (a) 
shows the full range of multiplicity, (b) shows expanded region of low multiplicity (ncAL < 20, nL0 < 10). 
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0) in 




(c) Oc^^- 



FIG. 2. The forward tile versus calorimeter tower multiplicity for (a) the central W boson sample, 
(b) the forward W boson sample, and (c) the Z boson sample. 

ment: < Et >= 35.2 versus 35.1 GeV, < >= 36.9 versus 36.5 GeV, and < Mr >= 70.4 versus 
72.5 GeV/c^, with uncertainties of about 3% on the latter values due to the hmited statistics for 
diffractive candidates (91 events). 

The fractions of W and Z boson events that contain forward rapidity gaps ("gap fraction") 
are extracted from fits to the data in Fig. 2. The non-diffractive (high multiplicity) background 
in the signal region is represented by a four-parameter polynomial surface, and the signal by a 
two-dimensional falling exponential as in Ref. [11]. Figure 4 shows the multiplicity distribution 
from Fig. 2(a), and the resulting fitted signal, fitted background, and the normalized distribution 
of pulls ( [data-fit] /-\/]V). Thex^/dof = 1.04 for this fit, and all other fits are of comparable quality. 

The fitting process is repeated over a systematically varied range of (lower limit of 2 tiles 
and upper limit ranging from 5 to 7 tiles) and ncAL (lower limit of 3 towers and upper limit ranging 
from 6 to 12 towers) to minimize any dependence of the signal on the chosen region. The result 
is a distribution of gap fractions, with the final signal defined by the mean of this distribution. 
The statistics in the Z boson sample are insufficient to perform an independent fit, so we use the 
background shape from the combined W boson sample scaled to the Z boson data to determine 
the diffractive Z boson signal. Varying the shape of the background samples used for the fit shows 
only small variations in the signal, well within the quoted uncertainties. 

To determine the final gap fractions, we correct the fitted values for residual contamination from 
multiple interaction events which were not rejected by the single interaction requirement. These 
events contribute only to the denominator of the gap fraction, resulting in a measured gap fraction 
that is lower than the correct value. This correction is determined using the full data sample with 
no single interaction requirement by comparing the predicted number of single interaction events 
(based on the instantaneous luminosity) with the observed number of events after the single 
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FIG. 3. Event characteristics for standard W boson events (left column, ncAL > 1), compared to 
diffractive W boson candidates (right column, ncAL = ?^l0 = 0)- The top plots compare electron Et, 
the middle plots show and the bottom plots compare the transverse mass (Mr) for the two cases. 




FIG. 4. The (a) ui^q versus ncAL distribution for the central electron W boson data from Fig. 2(a), 
and corresponding (b) fitted signal, (c) fitted background, and (d) normalized pull distributions. 
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interaction selection. This method demonstrates that our single interaction requirement is quite 
effective, and yields only an absolute correction of (0.09 ±0.05)% for the central electron W boson 
and the Z boson samples and a negligible correction for the forward electron W boson sample. 

Table 111 summarizes the final gap fractions obtained for the W and Z boson samples and their 
significances. The combined W boson sample has a gap fraction of 0.89 ±o;}7 % and a probability 
that the non-diffractive background fluctuated to the diffractive signal of 3 x 10~^^, corresponding 
to a significance of 7.5(T. The central W boson fraction (electron \ri\ < 1.1) of 1.08±Qjy% is greater 
than the forward fraction (1.5 < \ri\ < 2.5) of 0.64 ±o'ig %, unlike for jet events [11] (or typical 
diffractive expectations), which have a larger forward fraction. The Z boson sample has a gap 
fraction of 1.44 ±q52 %, with a significance of 4Aa. Uncertainties are dominated by those on the fit 
parameters. Additional small uncertainties from the dependence on the range of multiplicities used 
in the fits are added in quadrature. Potential sources of systematic error, such as the number of 
fit parameters, electron quality criteria, tower thresholds, and residual noise, yield only negligible 
variations in the gap fractions [18,21]. 



Sample 


Gap Fraction 


Significance 


Central W ^ eu {\rj\ < 1.1) 
Forward W ^ ev (1.5 < \ri\ < 2.5) 
Total W -^eu 
Total Z ^ e+e" 


(1.08 + 0.19 - 0.17)% 

(0.64 + 0.18-0.16)% 
(0.89 + 0.19-0.17)% 
(1.44 + 0.61 - 0.52)% 


1 X 10-^"^ {7.7a) 
6 X 10"^ (5.3cj) 
3 X 10"^^ (7.5(7) 
5 X 10"^ (4.4a) 



TABLE III. Measured gap fractions and probabilities for non-difFractive W and Z boson events to 
fluctuate and mimic diffractive W and Z boson production. 



We have thus far considered only non-diffractive W and Z boson events as the relevant back- 
ground to diffractive production. We now consider contamination from events other than the 
desired W ^ eu and Z ^ ee states, drawing primarily on the work of Ref. [20] . 

The largest background to W boson production is from multijet events in which one jet is 
misidentified as an electron, while another is measured incorrectly, thereby providing large Iprp. 
The fraction of fake W ^ ev events from multijet production was calculated [20] to be 0.046±0.014 
and 0.143 ± 0.043 of the total W ^ ev events for the central and forward electron samples, 
respectively. We use these fractions to determine the number of multijet events in our samples, 
and use measurements from Ref. [11] to obtain the number of diffractive events expected from this 
background. The results are shown in Table IV: for the central electron sample, a total of 0.88 
diffractive events are expected from the 401 multijet background events. Given that there are 8724 
events in the central electron W boson sample, with a measured diffractive fraction of 1.08 ±9^7%, 
we expect a total of 94±j^5 diffractive events. Recalculating the central W boson gap fraction after 
subtracting the multijet background gives a slightly higher value of 1.11% (93/8323). We note that 
this 3% change is an upper limit, because multijet background in events with single interactions 
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would likely be smaller than in the inclusive W boson sample due to smaller fluctuations expected 
in ^j,. The forward sample gives a negligible correction, since the gap fraction from diffractive W 
boson signal and multijet background are nearly identical. 



Sample 


Total 


Multijet 


Multijet 


Diffractive 


Diffractive 




Events 


Fraction 


Events 


Dijet Fraction 


Multijet Events 


central W 


8724 


0.046 ± 0.014 


401 


0.22 ± 0.05% 


0.88 


forward W 


3898 


0.113 ±0.013 


557 


0.65 ±0.01% 


3.6 



TABLE IV. The number of multijet background events in the diffractive W boson sample is calcu- 
lated. Then, given the number of multijet events in the sample and the diffractive dijet rate, the number 
of diffractive events expected from these background events is calculated. 



In addition to the multijet background, we consider background from misidentified Z boson 
events in which one electron is not detected. Again using methods from Ref. [20], we estimate 

94 ± 24 Z boson events in the combined W boson sample, with 1.35 of these being diffractive W 
boson candidates. Subtracting this background would result in a less than 1% correction, in the 
opposite direction from the multijet correction, since the diffractive Z boson signal is larger than 
that of the diffractive W boson. Finally, the background level from W ^ ru is expected to be 
small (about 2%), and we would expect the same gap fraction from W bosons that decay to r 
leptons as from the electron channel, therefore no correction is needed. 

Combining all these background sources yields a total background to diffractive W boson 
production of at most 2%, which is insignificant compared to the total 20% uncertainty, and we 
therefore do not apply any correction. Similar considerations for the diffractive Z boson sample 
yield at most a 4% background correction factor, which is again not significant, and consequently 
not applied. 

In this paper we have chosen to present the gap fraction, which is directly based on observable 
quantities. We thus avoid the reliance on potentially large model-dependent corrections. Therefore 
the measured gap fraction of 1.08 ±o % for central electron W boson events cannot be directly 
compared to the CDF measurement of (1.15 ± 0.55)% [16], which includes a correction factor 
derived from the POMPYT diffractive Monte Carlo [22] (based on the Ingelman-Schlein model) to 
attempt to account for how often a diffractive event does not yield a rapidity gap. They obtained 
a correction factor of 0.81 based on their Beam-Beam Counter (BBC) multiplicity, implying an 
uncorrected value of (0.93 ± 0.44)%, which is consistent with our measurement. However, this 
correction factor is quite different from that obtained by D0 and CDF using two-dimensional 
(luminosity counter and calorimeter multiplicity) methods subsequently adopted by both collab- 
orations to extract rapidity gap signals. We obtain a correction factor of 0.21 ± 0.04 from our 
diffractive W boson Monte Carlo using a quark structure for the pomeron, which compares well 
with the quark-structure-based correction for our central jet measurement (0.18) [18] and the CDF 
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correction factor for their diffractive 6-quark production (0.22) [7] and J/^^ production (0.29) [8]. 
Since there is no consensus on the correct model for describing diffractive data at the Tevatron, 
we feel that using POMPYT to correct the data is not advisable, but it should be noted that our 
corrected W boson gap fraction would be 5.1% according to this model, while there would be no 
correction needed for non-pomeron models such as SCI. 

Next, we calculate the ratio of the diffractive W and Z boson cross sections. In addition 
to intrinsic interest in this measurement, it is a potentially important input to the systematic 
uncertainty on the ratio R of the two cross sections [20]. We can write the diffractive cross section 
ratio Rd in terms of the gap fractions and the ratio R as follows: 

where we have substituted the measured gap fractions W^/W and ZiyjZ from this Letter, and the 
measured value R = 10.43 ± 0.15{stat) ± 0.20{syst) ± O.lO(iVLO) [20], which takes into account 
acceptance differences between the W and Z boson samples (assumed to be similar for diffractive 
and non-diffractive events). This value of Rd is consistent with the ratio for non-diffractive 
production. 

Finally, we measure the fractional momentum loss of the scattered proton ^ using the following 
equation [23]: 

e^^^E^T^e"^ (1-2) 

where Et^ and rji denote the transverse energy and pseudorapidity, respectively, of the observed 
particles. The rj of the outgoing scattered proton or antiproton (and the rapidity gap) is defined 
to be positive. As discussed in Ref. [11], Eq. 1.2 is particularly sensitive to particles emitted in 
the well-measured central region near the rapidity gap, while particles lost down the beam pipe 
at negative rj have negligible effect. Using a sample of POMPYT W boson events, where C, can 
be determined from the momentum of the scattered proton, we have verified that Eq. 1.2 is valid 
independent of pomeron structure. A scale factor 1.5 ± 0.3, derived by passing the Monte Carlo 
data through a full detector simulation, is used to convert ^ measured from all particles to that 
from just the electromagnetic energy depositions in the calorimeter [21]. Figure 5 shows the ^ 
distribution for the diffractive W boson candidate event sample with ncAL = "^10 = 0. The mean 
^ is 0.052 and most of the events have ^ < 0.1. Comparison of this ^ distribution obtained from 
calorimeter information with that from the measured proton using the upgraded D0 detector in 
Run II will give important insight into the nature of diffraction. 

In conclusion, we have observed diffractive W boson production with greater than 7a signifi- 
cance, shown that these diffractive W boson candidates have similar properties to non-diffractive 
ones, and measured the fraction of W boson events that are produced diffractively, in both the 
central and forward regions. We also have provided the first evidence for diffractive Z boson 
production. The extracted gap fractions have no model-dependent corrections, and are typically 
about 1%, far below expectations for a quark-dominated pomeron. We have obtained a ratio of 
diffractive W and Z boson cross sections consistent with the ratio for non-diffractive production. 
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FIG. 5. The ^ distribution for the combined central and forward electron W boson data with 
ncAh = ni^Q = 0, extracted from calorimeter information as described in the text. 



We have also measured the fractional momentum lost by the scattered proton and found that typ- 
ically less than 10% of the proton momentum takes part in the hard scattering, with an average 
of about 5%. 
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